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ABSTRACT 

The Rossi X-ray Timing Explorer has observed five outbursts from the transient 2.5 ms accretion- 
powered pulsar SAX J1808.4— 3658 during 1998-2008. We present a pulse timing study of the most 
recent outburst and compare it with the previous timing solutions. The spin frequency of the source 
continues to decrease at a rate of (— 5.5± 1.2) x 10~^^ Hz s~^, which is consistent with the previously 
determined spin derivative. The spin-down occurs mostly during quiescence, and it is most likely due 
to the magnetic dipole torque from a B = 1.5 x 10^ G dipolar field at the neutron star surface. We 
also find that the 2 hr binary orbital period is increasing at a rate of (3.80 ± 0.06) x 10~^^ s s~^, 
also consistent with previous measurements. It remains uncertain whether this orbital change reflects 
secular evolution or short-term variability. 

Subject headings: binaries: general — stars: individual (SAX J1808.4— 3658) — stars: neutron — 
stars: rotation — X-rays: binaries — X-rays: stars 



1. INTRODUCTION 

Ten years ago, observations of the recurrent X- 
ray transient SAX J1808. 4-3658 with the Rossi X- 
ray Timing Explorer (RXTE) revealed coherent 2.5 ms 
pulsations, establishing this source as the first 

known accretion-powere d millisecond pulsar (AMP; 

IWiinands fc van der KlisI [T998i) . Its immediate recogni- 
tion as the "missing link" in which old pulsars are spun 
up to millisecond periods by accretion in an X-ray bi- 
nary was only the beginning of what can be learned 
from this class of sources: among other areas of theo- 
retical interest, AMPs continue to reveal insights into 
the nature of magnetically channeled accretion flow, the 
torques that act upon neutron stars in low-mass X-ray 
binaries (LMXBs), and the orbital evolution of these sys- 
tems. This paper investigates the 2008 October X-ray 
outburst of SAX J1808.4— 3658, which touches upon all 
these questions. 

SAX J1808.4— 3658 was originally discovered by the 
BeppoSAX WFC monitor during a 1996 outburst 
(|in 't Zand et al.l 11998. 2001). It was the first LMXB 
known to exhibit all three modes of millisecond X- 
ray variability esta blished by the RXTE: accre tion- 
powered pulsations llWiinands fc van der Klisl [l99l. X- 
ray burst oscillations (Chakraba rtv et al.l2003l), andkilo - 
hertz quasi-periodic oscillations ( Wijnands et al.ll2003l ). 
The source has been extensively studied for over 1.7 Ms 
with the RXTE PGA during outbursts in 1998, 2000, 
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2002, 2005, and 2008, making it the best-studied AMP. 
Its outbursts rise from quiescence (Lboi 
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ICampana et~alll2002t iHeinke et ani2007f r to their peak 
fiuxes in «5 d, followed by w5 d at peak bolometric 
luminosities of iboi — (5-8) x 10"^^ erg s~^ (for a dis- 
ta nce of 3.5 kpc and the bolo metric correction given 
bv iGallowav fc Cummind l2006l ). ~15 d of decay back 
to near-quiescent fluxes, and a month or more of low- 
luminosity flaring (Hartjnan .ct al. 20 08ai . hereafter H08). 
These outbursts recur every 1.6-3.3 years. 

Orbital modulation of the persistent pulsa- 
tions establishes a 2.01 hr binary orbital period 
(jChakrabarty fc MorganI Il998l ). Its mass donor is 
likely an extremely low-mass (» 0.05 Mq) brown 
dwarf (|Bildsten fc Chakrabartv! 1200 1[ ). The system has 
been detected and extensively studied in the optical 



been detecte d and ext ensively stud i ed m the optical 
(iRoche et al l 17998: Gil es et all [l999l: iWang eiTall ^OOH : 
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rela- 
tively high optical luminosity during X-ray quiescence 
has led to speculation that the neutron star (NS) 
may be an active radio pulsar during these intervals 
(iHoiner et al.l 120011: [Blirderi et al.l 120031 : ICampana et all 
l2004f ). although rad io pulsations have not been de- 
tected (jBurgav et al.l [200 3) . Pulse timing of the four 
outbursts during 1998-2005 revealed that the spin 
frequency of SAX J1808.4— 3658 is decreasing during 
quiescence at a rate of = (—5.6 ± 2.0) x 10^^^ Hz s^^, 
while its orbital period is increasing at a rate of 
Po,h = (3.5 ± 0.2) X 10 
[2008h . 
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2. X-RAY OBSERVATIONS AND ANALYSIS 

Routine monitoring of the Galactic bulge region with 
the RXTE PCA revealed S AX J1808.4-3658 to aga in be 
in outburst on 2008 Sep 21 (jMarkwardt fc Swanlj |2008'). 
RXTE observed the source almost every day during the 
outburst, producing 270 ks of data with 122 /is time res- 
olution and 64 energy channels covering the 2-60 keV 
range of the PCA. Unfortunately, the total effective area 
of the PCA has decreased significantly since the earlier 
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outbursts: an average of only 1.8 of the 5 proportional 
counter units are on during these observations, resulting 
in a mean effective area of 2200 cm^ . 

Our analysis followed the techniques developed 
in H08 and applied to the earlier outbursts of 
SAX J1808.4-3658. We shifted the photon arrival times 
to the solar system barycenter using the optical posi- 
tion from H08, applied the RXTE fine clock correction, 
and filtered out data during Earth occultations and in- 
tervals of unstable pointing. We also removed data from 
1 min before to 10 min after the two thermonuclear X-ray 
bursts observed by the RXTE during this outburst. 

To measure the times of arrival (TO As) and fractional 
amplitudes of the persistent pulsations, we folded 512 s 
intervals of data using a preliminary timing model based 
on extrapolating forward the H08 results. The timing 
models were applied and fitted using the TEMPO pulse 
timing program, version 11.005,^ and assumed a circu- 
lar orbit and a fixed spin frequency. From the resulting 
folded pulse profiles, we measured the phases and am- 
plitudes of the fundamental and the second harmonic of 
the pulsation.* Higher harmonics had insufficient power 
to be useful for pulse phase timing. 

Special care must be taken when measuring the pulse 
TOAs for SAX J1808.4-3658. As in past outbursts, 
the harmonic phase residuals relative to our best timing 
models exhibit variability on time scales of ~10 hr and 
longer that is well in excess of the Poisson noise expected 
from counting statistics. As described in H08, we applied 
a frequency-domain filter that weighs the harmonics ac- 
cording to observed noise properties in order to gener- 
ate a minimum- variance estimator of the NS spin phase. 
For this outburst, the second harmonic had less excess 
timing noise, so it was weighted more heavily when esti- 
mating the overall pulse TOAs. With these pulse TOA 
estimates, we used TEMPO to fit an improved timing 
model. By iteratively re-folding the data, estimating the 
pulse TOAs, and fitting a new model, we converged on 
our best-fit timing model. 

Uncertainties in our frequency measurements were esti- 
mated using Monte Carlo simulations, following the pre- 
scription of H08. These simulations account for the im- 
pact of long-time-scale correlations in the pulse arrival 
times. To test for the presence of a frequency derivative, 
we fit one to the pulse TOAs and estimated its signifi- 
cance with our simulations. 

To measure the orbital period and its derivative, we 
connected the orbital phases of all the outbursts. The 
uncertainty in the 2008 orbital phase was estimated by 
measuring and accounting for the timing noise at scales 
of < Porb; as in previous outbursts (H08), the noise was 
close to white but somewhat higher than what would be 
expected form Poisson noise alone. 

3. RESULTS 

In most respects, the 2008 outburst of 
SAX J1808.4— 3658 continues the trends and pat- 
terns observed during the earlier outbursts. Figure [T] 
shows its light curve, pulse phases, and fractional 
amplitudes, along with six representative 2-15 keV pulse 
profiles. 

At |http: / /wwww.atnf .csiro. au / research /pulsar / tempo [ 
® Throughout this paper, we number the harmonics such that the 
fcth harmonic is k times the frequency of the 401 Hz fundamental. 
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Fig. 1. — The pulse profiles, light curve, phase residuals, and frac- 
tional amplitudes for the 2008 outburst of SAX J1808. 4-3658. The 
top panels show the 2-15 keV pulse profiles observed during the 
outburst. The profiles are background-subtracted and normalized 
such that the phase bins have a mean value of unity, and the error 
crosses in the lower right of each graph give the mean uncertainties. 
Each box shows two cycles, so each minor tick mark on the phase 
axis is 0.1 cycles. The phases are based on our best- fit constant- 
frequency spin model. The grey boxes in the light curve plot indi- 
cate the time intervals during which the profiles were observed. The 
second panel shows the light curve of the 2008 outburst, in black. 
The remarkably similar 2005 outburst is shown in grey for compar- 
ison. The stars indicate the times of thermonuclear bursts during 
2008. The third panel shows the pulse phase residuals relative 
to our best-fit constant-frequency model, v = 400.975210133 Hz. 
Blue points indicate the phase of the fundamental; red points in- 
dicate the second harmonic. Each point represents 512 s of timing 
data. The phase jumps after MJD 54760 led us to separate the 
phase measurements as earlier (circles) and later (squares), and 
the dashed and dotted lines show the respective phase models for 
each group. The bottom plot shows the fractional amplitudes of 
the fundamental and harmonic. Arrows give 90% upper limits. 
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The light curve of the main body of the outburst 
is remarkably similar the light curve observed in 2005, 
which is included in Figure [T] for comparison. The 
2008 outburst is slightly dimmer (maximum 2.5-25 keV 
flux of (1.74 ± 0.02) X 10"^ erg cm^^ g-^) and its 
light curve has two maxima, but the outburst rise, 
the length of the peak, and the slow decay down to 
Ri4 X 10"^*^ erg cm~^ s~^are nearly identical. 

The flaring tail of the new outburst is markedly dim- 
mer. During the two previous outbursts for which the 
tail was observed immediately following the main out- 
burst (2002 and 2005), the flares were quasiperiodic with 
a 4-6 d time scale, and the 2.5-25 keV flux varied be- 
tween ~10^^^ erg cm~^ s~^ and 4 x 10^^° erg cm~^ s~^. 
The source remained consistently detectable for at least 
a month following the main outburst. In contrast, the 
brightest observations during the tail of the 2008 out- 
burst only reach 2.5 x 10"^*^ erg cm~^ s~^, and during 
Oct 13-17, the source was not detectable at all with the 
RXTE (~10~^^ erg cm~^ s~^ sensitivity). This behavior 
is similar to the inte rmittent detections du ring the tail 
of the 2000 outburst (jWiinands et al.ll200lD . 

Most of the pulse profiles were similar to what we ob- 
served during the earlier outbursts: profiles 1-3 of Fig- 
ure [1] closely resemble profiles 1-3 in Figure 3 of H08. 
The evolution of the pulse shape progresses in the same 
way. As seen in the 2005 and (to a lesser extent) 2002 
outbursts, a symmetric pulse during the burst rise be- 
comes increasingly skewed during the burst peak. In the 
phase residuals, this pulse shape change manifests as a 
progressively increasing phase of the fundamental, while 
the phase of the second harmonic remains constant; this 
behavior can be seen during MJD 54732-54735 in the 
phase offset plot of Figure [TJ Since these pulse shape 
changes confuse measurement of the spin phase, we omit 
the outburst rise and early peak (points before MJD 
54735) from our timing models. This choice is admit- 
tedly somewhat arbitrary, as the pulse profile continues 
to change throughout the outburst, albeit more slowly. 
However, this same region was omitted from the analysis 
of the earlier outbursts in II08, so we remove it here so 
that our frequency measurements are directly compara- 
ble.9 

The clear double peaks of profile 4 are an exception 
to these similarities in the pulse shape evolution. This 
profile folds 8.2 ks of data during MJD 54743.2-54744.0. 
It occurs during the slow-decay stage of the outburst, at a 
2.5-25 keV flux of l.Ox 10~^ erg cm~^ s~^. In none of the 
earlier outbursts was a second peak so prominent, and 
the weak secondary peaks that were observed occured 
exclusively during the flaring tail, when the source was 
a factor of 2-3 dimmer. 

The most unusual timing feature of this outburst is 
the presence of abrupt shifts of the fundamental's phase 
during the outburst tail. These phase jumps are clear at 
MJD 54770 and 54774 in the phase offset plot of Figure[Tl 
at which times the fundamental phases are 0.2-0.3 cycles 
earlier than predicted by the constant-frequency model 
shown. The observations with earlier phases were folded 

^ Including the phases during the outburst rise decreases the 
frequencies cited throughout this paper by 20—30 nHz. Including 
this stage in the earlier outbursts causes their measured frequencies 
to shift similarly, so the value of the long-term spin down is not 
affected. 



and summed to produce profile 6, and their phases are 
marked by circles. The rest of the observations in the tail 
were folded to produce profile 5, and after MJD 54760 
their phases are marked by squares. While substantial 
phase shifts have been observed in the prev ious outbursts 
of SAX J1808.4-3658 (compare Fig. 2 of iBurderi et~an 
l2006l and Fig. 1 of H08), these jumps were particularly 
large and rapid. During the 23 hours between the obser- 
vations on MJD 54770 and 54771, the phase shifted by 
0.26 cycles; during the single observation on MJD 54775, 
the phase shifted by 0.20 cycles in 40 minutes. These 
shifts are due in part to differing pulse shapes. The cen- 
troid of the pulse — to which the fundamental is most 
sensitive — arrives appreciably earlier in profile 6, but 
the peaks of the two profiles are separated by a smaller 
difference. However, this nuance is lost unless many ob- 
servations are folded together (as in the displayed pro- 
files), and for most observations in the tail only the fun- 
damental is available for pulse timing. 

Given this limitation, the inclusion of these phase mea- 
surements of the fundamental in our timing models is 
subject to interpretation. If these points are all ex- 
cluded, the phase residuals for MJD 54735-54760 give 
a frequency of 400.97521012(3) Hz. (Parenthetical digits 
give 1 a estimates of the uncertainties.) If the circled 
points are used and the square points excluded, the tim- 
ing solution requires a large apparent spin up to account 
for the advancing phases in the outburst tail. Its fre- 
quency rises at a rate v = (1.1 ± 0.2) x 10^^^ Hz s~^ 
from a value of Vstavt = 400.97521002(5) at the start of 
the outburst. This fit is indicated by the dotted line in 
the phase offset panel of Figure[T] It does not completely 
eliminate the jumps: a substantial difference (0.16 cy- 
cles) remains between the points around MJD 54770 and 
54774. On the other hand, if we only include the square 
points in the outburst tail, the data are consistent with 
a constant spin frequency of 400.97521013(2) Hz. This 
selection eliminates all phase jumps greater than 0.1 cy- 
cles. Finally, if we include all the points in the tail, we get 
i> = 0.6(4) X 10-" Hz s-i and z^start = 400.97521009(6), 
although the phase jumps result in a poor fit and large 
uncertainties. 

Ultimately we are most interested in these frequency 
measurements to determine whether the long-term spin 
down observed in H08 has continued, and on this point 
the data are clear: regardless of our assumptions, the es- 
timated spin frequency at the start of the 2008 outburst 
is less than the 2005 frequency of 400.97521019(2) Hz 
(H08). Fitting a constant slope through the frequen- 
cies observed for the five outbursts yields a long-term 
spin-down rate of ;> = (—5.5 ± 1.2) x 10"^^ Hz s"^, as 
shown in the top panel of Figure O This rate and fig- 
ure use the 2008 timing model obtained when we only 
include the square points; we prefer this model on phys- 
ical grounds (see §4.ip . That said, using the f start fre- 
quency from the model for only the circled points does 
not greatly alter the result, slightly increasing its mag- 
nitude to -6.2 X IQ-^^ Hz s"^ However, the 2008 point 
becomes an outlier and the of the fit increases. 

These values are in excellent agreement with the spin 
down of (-5.6 ± 2.0) x lO'^^ Hz s'^ found for the 1998- 
2005 outbursts (H08). The frequency uncertainties esti- 
mated using our Monte Carlo simulations give = 9-7 
with 3 degrees of freedom. To account for this large x^. 
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TABLE 1 

Combined timing parameters for SAX J1808.4- 



3658 



, , I , , , I , , , I r 

V = (-5.5±1.2)xl0-i6 Hz s- 



Orbital period, Poi-b (s) ^ 

Orbital period derivative, Porb (10~^^ s s~^) 
Projected semimajor axis, Oxsini (light-ms) 
Time of ascending node, Tasc (MJD, TDB) 
Eccentricity, e (95% confidence upper limit) 
Spin frequency, u (Hz) ^ 

Spin frequency derivative, u (10~^® Hz s~^) 



7249.156980(4) 
3.80(6) 
62.812(2) 
52499.9602472(9) 
< 1.2 X 10"* 
400.975210240(9) 
-5.5(1.2) 



Porb and u are specified for the time Tasc (2002 Aug 14.0). 

we scaled up the frequency uncertainties such that the re- 
duced is unity when calculating the uncertainty for i/. 
The fit is not particularly sensitive to the source position: 
shifting it by 1 (7 along the ecliptic (using the uncertainty 
from the optical position in H08) changes the spin down 
by 0.2 X 10~^^ Hz s~^. No position significantly decreases 
the of the linear trendline. 

The farthest outlier from a constant frequency deriva- 
tive is the 2000 outburst, of which we only observed the 
flaring tail. As noted in H08, this frequency measure- 
ment may be complicated by small systematic differences 
in the observed pulse frequency between the main out- 
burst and its tail, which were observed during the 1998 
and 2002 outbursts. Removing this outburst lowers the 
to 1.9 (2 degrees of freedom), consistent with a con- 
stant derivative, and changes the long-term spin down to 
(-6.5 ±0.8) X 10-1^ Hz s"i. 

The measured time of ascending node^*^ continues 
to advance relative to a constant orbital period. A 
quadratic provides a good fit (x^ = 2.9 with 2 degrees 
of freedom), consistent with a constant orbital period 
derivative of Porb = (3.80 ± 0.06) x lO'^^ s s'^. This 
measurement is 1.5 a great er than the va lue in H08. It 
is in good agreement with iBurderi et al.l (|2009l) , which 
measured the orbital phase of the 2008 outburst using 
XMM-Newton data. All other orbital parameters differ 
by less than 1 a from the values listed in H08. Table [T] 
summarizes all the parameters for the long-term pulse 
timing of SAX J1808.4-3658. 

Proper interpretation of the phase residuals hinges on 
an accurate handling of the orbital parameters, particu- 
larly the orbital period, due to the sparsity of data in the 
outburst's tail. From the orbital parameters derived by 
connecting the advancing orbital phases of the 1998-2005 
outbursts (H08), we extrapolate a 2008 orbital period of 
7249.15763(6) s. Including the 2008 outburst's phase, 
this method gives 7249.15771(2) s. Fitting the orbital 
parameters using the fundamental for the first half of the 
ouburst only (cf. iPatruno et al.l [2008( 1 or using only the 
second harmonic give compatible periods of 7249.160(3) s 
and 7249.153(5) s, respectively. However, fitting to the 
phase of the fundamental over the entire outburst yields 
a period of 7249.174(2) s, 9 cr higher than the far more ac- 
curate values derived through phase connection. Because 
the phase measurements in the outburst tail are sparse 
due to the source's highly variable flux, their coverage 
of the orbital phase is limited, and TEMPO attempts to 
"correct" for the observed jumps by over- fitting Porb- 



The ascending node is the point of zero degrees true longitude, 
which is the same as zero degrees mean longitude for a circular orbit 
like this one. 
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Fig. 2. — The changing spin and orbital timing of 
SAX J1808.4-3658 during the 1998-2008 outbursts. Top: 
constant-frequency measurements of each outburst, showing the 
spin down of the star. The frequencies are relative to ug = 
400.97521000 Hz. The x's mark what the frequencies would be if 
the fit source position differed from the actual position by 2 cr along 
the ecliptic plane in the direction of increasing RA. Bottom: the 
observed times of ascending node, relative to the expected times for 
a constant period. The Tasc of each outburst comes progressively 
later, indicating a positive orbital period derivative. 



4. DISCUSSION 

4.1. The double-peaked pulse profiles and phase jumps 

The presence of double-peaked pulse profiles during 
the 2008 outburst was notable both because two peaks 
were so well defined (compare with the weaker secondary 
peaks shown in Fig. 3 of H08) and because they occurred 
while the outburst was still relatively bright: 2-3 times 
more luminous than when double-peaked profiles were 
first seen during the 2002 and 2005 outbursts. H08 spec- 
ulated that the appearance of a second peak may be 
caused by the recession of the accretion disk as the accre- 
tion rate drops, reveali n g the star's antipodal hot spot. 
Ilbragimov fc PoutanenI (|2009( ) went further, estimating 
from magnetic accretion theory and geometric consider- 
ations that the appearance of the antipodal hot spot at 
fluxes of Pboi ~ 0.8 X 10~^ erg cm"^ s^^ during 2002 and 
2005 would imply a surface magnetic field of B w (0.4- 
1.2) X 10^ G. 

The detection of double-peaked profiles at Pboi = 2.1 x 
10~^ erg cm~^ s~^ during 2008 requires this field esti- 
mate to be increased somewhat, although its weak depen- 
dence on flux (oc F^/^ in eq. 14 of llbragimov fc PoutanenI 
[20091) means that is still not in conflict with the B < 
1.5 X 10® G limit implied by the long-term spin down 
(H08; see also i i4.2p . More troublesome is the question 
of why the second spot does not remain visible as the 
accretion rate continues to drop. It is difficult to imag- 
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ine a change in the geometry or inner radius of the disk 
that temporarily reveals the antipodal hot spot but has 
no impact at all on the accretion rate and flux, as is the 
case here. More exotic explanations, such as a tempo- 
rary bifurcation of the accretion column to produce two 
hot spots, may be needed; simulations have suggested 
that such configurations are easil y possible if high er- 
order field components are present ()Long et al.ll2008D . 

The large phase jumps in the tail of the 2008 outburst 
are the most unusual feature of the pulse profiles, with 
the phase of the fundamental in one case jumping by 
0.20 cycles in only 40 minutes. Clearly these are not 
changes in the rotational phase of the NS, which compli- 
cates their inclusion while constructing a timing solution. 
We considered two extremes: that the earlier points re- 
flected the NS spin phase, which necessitated a spin up of 
1.1(2) X lO^^'^ Hz s~^; and that the later points tracked 
the spin, in which case a constant-frequency timing so- 
lution provides a good fit. 

We prefer the latter case on physical grounds. 
From standard magnetic accretion torque theory (e.g., 
iPringle fc Reesl[l972r ). an accretion rate of M can impart 

a maximum torque of 27r/;> < M(G'Afrco)^^^, where r^o 
is the co-rotation radius of the NS and / is its moment 
of inertia. Estimating M = L\^o\{GM/ R)^^ , we find the 
maximum spin up of 

z><6 X 10-14 f -"--v ) fB^\ X 

— \ 1x10 -J erg cm^ s ^ I \2.12 / 

(3.5 kpc) ( 10^^ g cm^ ) (10 km) ^ 

(o^)"''(iofm)"'' Hzs- (1) 

Here we use the bolometric correct ion Cboi and distance d 
from lGallowav fc Cummin^ (|2006f ) and canonical values 
for the NS parameters. Integrating this expression for 
the observed flux over the 40 d outburst outburst, we 
can approximate a physical upper limit for the total spin 
up of Ai' < 0.12 /iHz. For the earlier points to reflect the 
spin phase, the frequency shift would need to be Av = 
1.1 X 10"" Hz s-i ■ 40 d = 0.38 /iHz. Unless the source's 
luminosity underestimates the accretion rate by a factor 
of a few, these points cannot reflect the NS spin phase. 

Of course, it is possible that the NS spin phase lies 
somewhere between these two extremes. Based on the 
data available, we cannot exclude the possibility of a 
spin up during the outburst. Indeed, though the con- 
fidence intervals in 1998-2005 include zero, they are bi- 
ased toward a positive O (H08), suggestive that a small 
sytematic spin up may be present. 

4.2. Interpreting the spin frequency derivative 

The observed long-term spin down of 
SAX J1808.4— 3658 produces a loss of rotational 
energy at a rate of Esd — ^tt^Ivv — —9 x 10"^^ erg s"^, 
assuming a canonical value of / = lO'*^ g cm^ for the NS 
moment of inertia. This spin down acts on the NS during 
quiescence (H08), and it is consistent with a constant 
rate over the ten years of observation. We consider three 
possible torques: magnetic dipole radiation, the cen- 
trifugal expulsion of matter from the inner accretion disk 
by the magnetic field (the so-called "propeller effect"), 
and gravitational radiation due to a mass quadrupole 



moment of the NS. We assume that these independent 
mechanisms contribute additively to the observed 
spin-down luminosity, E^d = ^^dipoio + £^prop + i^gr- 

The presence of pulsations during the outburst peaks 
requires B > 0.4 x 10^ to truncate the accretion 
flow above the NS surface (jPsaltis fc Chakrabartvl[l999t 
H08). This lower limit implies a magnetic dipole contri- 
bution of ii^dipoic/i^sd > 10%. If the field strength is only 
a few times higher, B — l.bx 10* G, then magnetic dipole 
radiation would account for the entirety of the long-term 
spin down. These tight limits on B leave a very small 
window for other mechanisms. Their contribution would 
have to be within a single order of magnitude of the con- 
tribution from magnetic dipole radiation. Given the wide 
physical ranges available to these independent spin-down 
torques, such a coincidence is unlikely. Additionally, the 
long-term spin down is consistent with a constant rate, 
as would be expected if i^dipoie dominates; the constancy 
of iJprop and is very model dependent. 

Independent estimates of the magnetic field strength 
suggest t hat B is toward the u pper end of our allowed 
interval. iGilfanov et al.l (|1998f ) noted that if the rapid 
dimming at the end of the main outburst is due to the 
onset of the propeller effect (as suggested by other out- 
burst properties, e.g., iHartman et al.l [2b08b[ ) . it would 
imply i? ~ 1 x 10* G. As noted in the previous section, 
the appearance of double-peaked pulse profiles in the tail 
sugge st a similar field strength (jlbragimov fc PoutanenI 
l2009f ). The detection of a relativistically broadened iron 
line during the 2008 outburst also suggest s a higher field 
of (3 lb 1) x 10* G (iCackett et al.l 12001 iPapitto et all 
I2009f) . For all these reasons, it is likely that magnetic 
dipole radiation provides the dominant torque. 

4.3. Interpreting the orbital period derivative 

The cause of the orbital period deriva, tive remains 
a subject of debate. H08, iDi Salvo et all ([2008,), and 
iBurderi et al.l |2^009) show that it cannot be accounted 
for solely by conservative transfer from the low-mass 
companion to the NS. We consider two alternatives: the 
ejection of mass from the system due to the ablation of 
the companion, and the influence of short-term effects 
that can exchange angular momentum between the com- 
panion and the orbit. 

Ejection of matter due to the ablation of the compan- 
ion. If magnetic dipole torque is the predominant con- 
tributor to the spin down of SAX J1808.4— 3658, as sug- 
gested in the previous section, then the source may be- 
have like a rotation-powered pulsar during quiescence, 
producing radio pulsations and a particle wind. While 
no ra dio pulsations have been detected (jBurgay et al.l 
[200l . the upper limits are not particularly constrain- 
ing. However, during quiescence the syste m is substan- 
tially brighter in the optical than expected (jHomer et al.l 
I2001h . This optical excess may be due to the heating of 
the companion by a particle wind (|Burderi et al.|[2003l ). 

If this particle wind ejects a large amount of matter 
from the system, it could explain the rapid orbital pe- 
riod change. For a companion modeled with a polytropic 
index of n = 3/2 (representative of a degenerate or fully 
convective star), a mass loss of ~10~^ Mq yv~^ would be 
needed to yield the observed Porb (|Di Salvo et al.| [2008l). 

Is the spin-down luminosity of E^d ~ —9 x 10'^^ erg s""'^ 
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sufficient to drive such mass loss? Assuming that it 
is radiated isotropicaUy, the incident luminosity onto 
the companion is i?abi = —j{Rc/a)^Esd, where Rc is 
the companion radius and a is the binary separation. 
Since the companion fills its Roche lo be, Rda » 0.1 5 
for a mass ratio of g « 0.05/1.4 (jEggletonl I1983D . 
The resulting luminosity onto the companion is E'abi = 
2.1 X 10'^^ erg s~^. Assuming perfect efficiency, that 
incident luminosity can drive a mass loss of Mc = 
-E^uRc/GMc = -1.2 X 10"^ Mq yr-i. Thus this sce- 
nario for orbital evolution is energetically feasible. 

Nevertheless, some serious problems with this model 
remain. A particle wind sufficient to eject over 95% of 
Mc from the system, as necessitate d by this model, would 
also drive away the accretion disk. iDi Salvo et al.l (|2008l ) 
suggest that outbursts therefore occur when the pres- 
sure of matter driven from the companion temporarily 
exceeds the radiation pressure. However, the outburst of 
SAX J1808.4-3658 and the other AMPs look quite sim- 
ilar to some of the outbursts form other LMXBs, both 
in time scales and shape. The disk instability model has 
been generally successful in explaining t he light curves 
of so ft X-ray transients in gen eral (e.g.. | King fc Ritted 
|1998[) and AMPs in particular ijPowell et al.ll2007D . It is 
difficult to reconcile these similarities in light curve shape 
and time scale with a very different mechanism governing 
the outbursts of SAX J1808.4-3658. 

Short-term changes in the orbital period. Alternatively, 
the observed Porb may not be representative of the sec- 
ular evolution, but is instead due to some short-term 
interchange of angular momentum between the compan- 
ion and the orbit. Similarly large orbital period deriva- 
tives are present in two so-called "black widow" mil- 
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lisecond radio pulsars, and these period derivatives have 
been observed to change sign on a ~10 yr time scale 
(jArzoumanian et al.l [l994l : iDoroshenko et all 1200 1[ ). It 
has been suggested that tidal dissipation and magnetic 
activity in the companion are responsible for the or- 
bital variability, requiring that the companion is at least 
partiall y non-degenerate, convec t ive, and magnetically 
active (* Arzoumani an et al.l Il994t lApplegate fc ShahamI 
119941 : [Doroshenko et all 120011 ) Eclipse timing of the 
LMXB EXO 0748-676 also shows rapid variation in the 
orbital period; explanations for the exchange of angular 
momentum include a "wobble" of the companion's iner- 
tial moments due to differential rotation and conv ection 
produced by its uneven heating (iWolff et al.ll200l) . 

However, in all the above systems, the orbital period 
derivatives change over <10 yr. While this does not 
definitively prove that the Porb of SAX J1808.4-3658 is 
secular, it would make this source an exception. In con- 
trast, the spin-down rate of the ablation model should 
be constant, so it is supported by the data from the 2008 
outburst. 



We are grateful to Jean Swank and the RXTE op- 
erations team at NASA Goddard Space Flight Center 
for their help in scheduling these observations. We 
also thank Tiziana Di Salvo, Luciano Burderi, Dun- 
can Galloway, and Mike Wolff for useful discussions. 
We thank the referee for useful suggestions. This work 
was supported in part by NASA grants NNX07AP93G 
and NNX08AJ43G, awarded to MIT through the RXTE 
Guest Observer Program and the Astrophysics Data Pro- 
gram. 



Giles, A. B., Hill, K. M., & GreenhiU, J. G. 1999, MNRAS, 304, 47 
Gilfanov, M., Revnivtsev, M., Sunyaev, R., & Churazov, E. 1998, 
A&A, 338, L83 

Hartman, J. M., Patruno, A., Chakrabarty, D., Kaplan, D. L., 

Markwardt, C. B., Morgan, E. H., Ray, P. S., van der KUs, M., 

& Wijnands, R. 2008a, ApJ, 675, 1468 (H08) 
Hartman, J. M., Watts, A. L., & Chakrabarty, D. 2008b, ApJ, 

submitted, arXiv: 0809.3722 
Heinke, C. O., Jonker, P. G., Wijnands, R., & Taam, R. E. 2007, 

ApJ, 660, 1424 

Homer, L., Charles, P. A., Chakrabarty, D., & van Zyl, L. 2001, 

MNRAS, 325, 1471 
Ibragimov, A. & Poutanen, J. 2009, MNRAS, submitted, arXiv: 

0901.0073 

in 't Zand, J. J. M., Cornelisse, R., Kuulkers, E., Heise, J., Kuiper, 

L., Bazzano, A., Cocchi, M., MuUer, J. M., Natalucci, L., Smith, 

M. J. S., & Ubertini, P. 2001, A&A, 372, 916 
in 't Zand, J. J. M., Heise, J., MuUer, J. M., Bazzano, A., Cocchi, 

M., Natalucci, L., & Ubertini, P. 1998, A&A, 331, L25 
King, A. R. & Ritter, H. 1998, MNRAS, 293, L42 
Long, M., Romanova, M. M., & Lovelace, R. V. E. 2008, MNRAS, 

386, 1274 

Markwardt, C. B. & Swank, J. H. 2008, The Astronomer's 
Telegram, 1728 

Papitto, A., di Salvo, T., D'Ai, A., laria, R., Burderi, L., Riggio, 

A., Menna, M. T., & Robba, N. R. 2009, A&A, 493, L39 
Patruno, A., Hartman, J. M., Wijnands, R., van der Klis, M., 

Chakrabarty, D., Morgan, E. H., & Markwardt, C. B. 2008, The 

Astronomer's Telegram, 1760 
Powell, C, Haswell, C., & Falanga, M. 2007, MNRAS, 374, 466 
Pringle, J. E. & Rees, M. J. 1972, A&A, 21, 1 
Psaltis, D. & Chakrabarty, D. 1999, ApJ, 521, 332 
Roche, P., Chakrabarty, D., Morales-Rueda, L., Hynes, R., Slivan, 

S. M., Simpson, C, & Hewett, P. 1998, lAU Circ, 6885 



A Decade of Timing SAX J1808.4-3658 



7 



Wang, Z., Bassa, C, Gumming, A., & Kaspi, V. M. 2008, ApJ, 

accepted, arXiv: 0812.2815 
Wang, Z., Chakrabarty, D., Roche, P., Charles, P. A., Kuulkers, 

E., Shahbaz, T., Simpson, C, Forbes, D. A., & Helsdon, S. F. 

2001, ApJ, 563, L61 
Wijnands, R., Mendez, M., Markwardt, C, van der Klis, M., 

Chakrabarty, D., & Morgan, E. 2001, ApJ, 560, 892 



Wijnands, R. & van dcr Klis, M. 1998, Nature, 394, 344 
Wijnands, R., van der Klis, M., Homan, J., Chakrabarty, D., 

Markwardt, C. B., & Morgan, E. H. 2003, Nature, 424, 44 
Wolff, M. T., Hertz, P., Wood, K. S., Ray, P. S., & Bandyopadhyay, 

R. M. 2002, ApJ, 575, 384 



